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Abstract. We investigate theoretically and experimentally how the hydrody- 
namically correlated lateral motion of particles in a suspension confined between 
two surfaces is affected by the suspension concentration. Despite the long range of 
the correlations (decaying as 1/r 2 with the inter-particle distance r), the concen- 
tration effect is present only at short inter-particle distances for which the static 
pair-correlation is nonuniform. This is in sharp contrast with the effect of hydro- 
dynamic screening in unconfined suspensions, where increasing the concentration 
changes the prefactor of the large-distance correlation. 



1. Introduction 

The Brownian motion of colloid particles and macromolecules is correlated through 
hydrodynamic interactions, i.e., flows that the motion induces in the host liquid ^12- 
In an unconfined suspension [3] these correlations decay with inter-particle distance 
r as 1/r. They are positive, i.e., particles drag one another in the same direction. 
The long range of the interaction leads to strong many-body effects, manifest in an 
appreciable dependence of transport coefficients on concentration. A concentration 
effect of particular interest is hydrodynamic screening [2] , which sets in over distances 
much larger than the inter-particle distance and renormalizes the prefactor of the 
~ 1/r pair interaction. 

Colloids may sometimes be spatially confined by rigid boundaries, as in porous 
media, biological constrictions, or microfluidic devices. Attention has been turned 
recently toward the dynamics of such confined suspensions 0J-^2] due to the 
emergence of microfluidic applications and the development of new visualization and 
manipulation techniques (digital video microscopy |13j and optical tweezers |14p. 
allowing to study dynamics at the single- and few-particle level. The dynamics of 
confined suspensions have been studied also by computer simulations |15j-|17) and 
have renewed the interest in related hydrodynamic problems |18l 

We have recently demonstrated the dramatic effect that confinement between two 
flat surfaces has on the hydrodynamic pair interaction at large inter-particle distances 
|H3 El- Such a quasi- two-dimensional (Q2D) suspension is portrayed in figure ^ 
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Despite the confinement, the hydrodynamic interaction is still long-range, decaying 
as 1/r 2 rather than 1/r. There is "anti-drag" in the transverse direction, i.e., the 
correlation between two particles moving perpendicular to their connecting line is 
negative. Arguably the most striking finding, however, is that the concentration 
of the suspension has no effect on the large-distance correlation, i.e., there is no 
hydrodynamic screening. These three properties are in stark contrast with their 
unconfined counterparts mentioned above. 

In the current paper we review these results and then extend the analysis to 
shorter inter-particle distances where the static pair-correlation of the suspension 
plays a crucial role. Throughout the paper we present the theoretical analysis along 
with the corresponding experimental results. Section [2] introduces the model system 
and the corresponding terminology, and section describes the experimental setup 
and methods. In section 0] we briefly discuss the dynamics of a single particle and 
in section [S] address the hydrodynamic interaction between two isolated particles. 
Sectional which constitutes the main part of the current work, deals with the three- 
body correction to the pair-interaction at finite concentration. Finally, in section|7|we 
conclude and discuss the results. 






Figure 1. (a) Optical microscope image of the experimental Q2D suspension at 
area fraction <j> = 0.338. (b) Schematic view of the system and its parameters. 



2. Model system 

The geometry considered in this work is depicted in figure ^b). Identical, spherical 
particles of radius a are suspended in a liquid of viscosity rj and temperature T, 
confined in a slab of width w between two planar solid surfaces. The x and y axes are 
taken parallel, and the z axis perpendicular, to the surfaces, where z = is the mid- 
plane. The particles behave as hard spheres with no additional equilibrium interaction. 
For simplicity we consider cases where particle motion is restricted to two dimensions, 
i.e., to a monolayer lying at the mid-plane. We use the notation r(p, z) for three- 
dimensional position vectors, where p{x, y) is a two-dimensional position vector in the 
monolayer. The area fraction occupied by particles is denoted by 4>. The Reynolds 
number is very low, and the hydrodynamics, therefore, are well described by viscous 
Stokes flows [3] . The confining boundaries are impermeable and rigid, imposing no-slip 
boundary conditions on the flow. 

We characterize the dynamic pair correlation between two particles by the 
coupling mobilities £>£ T (p) as functions of the inter-particle distance p. These are 
the off-diagonal terms in the mobility tensor of a particle pair, i.e., the proportionality 
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coefficients relating the force acting on one particle with the change in velocity of the 
other. The two independent coefficients, -B£ and B^, correspond, respectively, to the 
coupling along and transverse to the line connecting the pair. The dynamic correlation 
between two Brownian particles is similarly characterized by two coupling diffusion 
coefficients, _D£ T (p), which, due to the Einstein relation, are simply related to the 
coupling mobility coefficients via the thermal energy, £>£ T (p) = k^,TB^ T (p). [In [TH] 
four coefficients were considered, T , whose relation with the ones considered here 
is Dt y = (D+ T - Dl T )/2] 

3. Experimental setup 

The experimental system consists of an aqueous suspension of monodisperse silica 
spheres (diameter 2a = 1.58±0.04 pm, density 2.2 g/cm 3 , Duke Scientific), undergoing 
Brownian motion while being tightly confined between two parallel glass plates in a 
sealed thin cell (figure P). The inter-plate separation is w = 1.76 ± 0.05 pm, i.e., 
slightly larger than the sphere diameter, 2a/w ~ 0.90. Digital video microscopy and 
subsequent data analysis are used to locate the centres of the spheres in the field of 
view and then extract time-dependent two-dimensional trajectories. Details of the 
setup and measurement methods can be found elsewhere [20] ■ Measurements were 
made at four values of area fraction, <p = 0.254, 0.338, 0.547, 0.619 ± 0.001. (Larger 
area fractions could not be checked because the suspension began to crystallize \21\.) 
From equilibrium studies of this system |21| we infer that, for the purpose of this 
study, the particles can be regarded as hard spheres. 

The coupling diffusion coefficients as functions of the inter-particle distance are 
directly measured from the tracked trajectories as 

Di{p) = (x 1 (t)x 2 (t)) p /(2t), D c T (p) = ( yi {t)y 2 {t)) p /{2t), (1) 

where Xi(t) and yi(t) are the displacements of particle i of the pair during a time 
interval t along and transverse to their connecting line, respectively. The average 
(}p is taken over all pairs whose mutual distance falls in a narrow range (±0.09 /im) 
around p. 

4. Single particle 

Consider a single particle whose centre, lying on the mid-plane between the two 
confining surfaces, is defined as the origin. A force fi is applied to the particle in 
the i direction parallel to the surfaces (i — x,y). As a result, the particle moves with 
velocity 

u\i = B s f u = B [{a/w)A s (a/w)}fn, (2) 

where B s is the self- mobility of the particle in the given geometry and Bq = (fiir-na^ 1 
is its self-mobility in an unconfined liquid. Alternatively, we may consider a free 
Brownian particle. Its mean-square displacement during a time interval t will be 

(p\t)) = AD [(a/w)A s (a/w)]t, (3) 

where D Q — k^TB^. The dimensionless factor (a/w)A s (a/w), representing the effect 
of confinement, becomes unity in the limit a/w <C 1. Approximate expressions for 
this factor for larger values of the confinement ratio a/w were the subject of many 
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previous works (see [3] and references therein), and their validity has been confirmed 
in recent experiments [H]-|E]- 

The particle motion makes the surrounding liquid flow. At distances much larger 
than w the flow velocity is given by ^3 E| 

Vi(r) = (a/w)B A ij (r)f lj 

A l3 (p, z)= - A ^ Ui - H{z/w), (4) 

P V P J 

where i,j = x,y, the vertical component v z is exponentially small in p/w, and Ao is 
a dimensionless prefactor of order f dependent only on the confinement ratio a/w. 
The flow field JIJ can be viewed as produced by an effective two-dimensional mass 
dipole (source doublet), located at the particle centre and oriented in the j direction 
|f 8U22llTTJllll| . The flow has a perpendicular profile H(z/w) which vanishes on the two 
confining surfaces, H{— 1/2) = H(l/2) = 0, and is normalized to f at the mid-plane, 
H(0) = 1. 

As argued in ^] ^] > t ne f &r now nas t ne dipolar shape |@} regardless of particle 
size. Changing the confinement ratio a/w merely modifies the effective mass-dipole 
strength, i.e., the prefactor Xo(a/w). In the limit a/w — » one finds Ao = 9/16 |18|. 
As will be shown below, we find for a/w ~ 0.45 Ao ~ 0.49. Since the maximum 
possible confinement ratio is a/w = 1/2, it is inferred that the function Xo(a/w) 
actually changes very little with a/w. This weak dependence on the confinement ratio 
can be understood in terms of mutual compensation of two opposing effects: when 
a/w is larger, on one hand, the particle displaces a larger liquid volume as it moves, 
yet, on the other hand, its self-mobility decreases [i.e., (a/w)A s (a/w) gets smaller]. 

5. Pair interaction 

Let us introduce a second particle whose centre lies at r = (p, 0). The particle is 
torque-free. It is also force-free in the xy plane. Being confined to the mid-plane, it 
evidently cannot be assumed force-free in the z direction. However, the flow J3J does 
not exert any force in that direction on a mid-plane-placed particle. Hence, as long as 
we restrict the discussion to such flows, the particle can be regarded as force-free. It 
is thus entrained by the flow with velocity |3] 

U2l = ^a?J A dr "^ r ") =( 1 ~ |J) °)' ^ 

where the integration is over the particle surface A. In obtaining the last equality 
we have assumed a parabolic (Poiseuille) vertical profile, H(z/w) = 1 — 4z 2 /w 2 |18|. 
Equation JSJ is, in fact, a manifestation of Faxen's first law |3] as applied to the 
Q2D geometry. For a force-free particle, Faxen's law yields u = v + (a 2 /6)V 2 v. The 
dipolar flow 10} has (d xx + d yy )v — and <9 zz | z =o v = [w~ 2 H"(0)/H(0)]v, from which 
equation (JSJ readily follows. 

Substituting equation Q in equation JSJ, we obtain a relation between the force 
fi exerted on the first particle and the resulting velocity change 112 of the second 
particle, 

u 2l = [1 - Aa 2 /(3w 2 )](a/w)B A l3 {p,0)f l3 , (6) 

thus identifying the off-diagonal coupling terms of the pair- mobility tensor as (p) = 
[1 — Aa 2 /(3w 2 )](a/w)B Aij(p,0). To get the coupling mobility along the line 
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connecting the pair, we consider the relation between fx x and u 2x for p — px. Using 
equation (QJ we get 

Bl{p) = [1 - ±a 2 /(3w 2 )](a/w)B Q A xx (px 1 0) = (a/w)B A L (p) 
A L (p) = Xw 2 /p 2 . (7) 

Similarly, for p = py the relation between f\ x and uix yields the coupling mobility 
transverse to the connecting line, 

S|(p) = [1 -4a 2 /(3w 2 )}{a/w)B A xx (py,0) = {a/w)B A T (p) 

A T (p) = -Xw 2 /p 2 . (8) 

In equations Q and © we have defined 

X(a/w) = [1 - 4a 2 /(3u> 2 )]A (a/w). (9) 

This is a refinement of our previous analysis llOHllj . In a Q2D geometry the so-called 
stokeslet approximation, equating the coupling mobility with the flow velocity per unit 
force, strictly holds only when a is much smaller than both the inter-particle distance 
p and the slab width w, whereupon A ~ A — 9/16. If a is much smaller than p but 
comparable to w then, no matter how large the inter-particle distance may be, we have 
A < Aq. This is due to the second term in Faxen's law, ~ a 2 V 2 v. In an unconfined 
liquid it contributes a negligible correction to the stokeslet limit, 0(a 2 /r 2 ), whereas 
in Q2D the newly introduced length w leads to a significant correction of 0(a 2 /w 2 ). 

As found in equations 10 and (JHJ), the pair hydrodynamic interaction in Q2D 
is very different from its unconfined counterpart. The decay with distance is faster, 
~ 1/r 2 instead of ~ 1/r, yet the interaction is still long-range |23l I22j . (Its decay, 
in fact, is slower than near a single surface, where the interaction falls off as 1/r 3 
The transverse coupling is negative, i.e., particles exert "anti-drag" on one 
another as they move perpendicular to their connecting line. (In the unconfined case 
one has A T = A L /2, both coefficients being positive.) These properties are direct 
consequences of the far flow field |@}. The 1/p 2 decay is that of the flow due to a 2D 
mass dipole. The negative transverse coupling is a result of circulation flows in the 
dipolar field [TO]. 

For a pair of force-free Brownian particles, the hydrodynamic coupling derived 
above implies that their motions be correlated according to equation 0, with 
D| T (p) = (a/w)£>oAL,T(p). Thus, the dimensionless couplings Al,t(p) can be 
directly measured from the statistics of particle trajectories. These measurements for 
different area fractions 4> are presented in figure [21 The large-distance behaviour for 
all 4> values fits well the ±\w 2 /p 2 dependence of equations and (|SJ) with A ~ 0.36. 
The negative sign of Ax confirms the predicted "anti-drag" between particles located 
transverse to the direction of motion. Substituting A ~ 0.36 and a/w ~ 0.45 in 
equation we find A (0.45) ~ 0.49, which is only slightly smaller than the value for 
a vanishing confinement ratio, Ao(0) = 9/16 ~ 0.56. 

6. Concentration effect 

As the area fraction cf> is increased, the pair hydrodynamic interaction should become 
affected by the presence of other particles. In unconfined suspensions hydrodynamic 
screening sets in at distances much larger than the typical inter-particle distance and 
renormalizes the interaction by a concentration-dependent prefactor 2 . Similarly, one 
expects the pair interaction in Q2D to have the form of equations J7J) and © yet with 
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Figure 2. Longitudinal (Al, circles) and transverse (At, squares) coupling 
diffusion coefficients as a function of inter-particle distance p. The coefficients are 
scaled by Dga/w and the distance by w. Area fractions are <j> = 0.254 (a, red), 
0.338 (b, orange), 0.547 (c, green), and 0.619 (d, blue). All data are redrawn 
in (e), demonstrating a concentration-independent collapse of the large-distance 
measurements onto the same two curves. Dashed lines are a fit to ±\/(p/w) 2 
with the same value of A = 0.36 for all panels. 



a modified, ^-dependent prefactor. We find, however, that this is not the case, as is 
clearly demonstrated in figure EJe). 
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Consider again particle 1, located at the origin and exerting a force fi parallel 
to the confining surfaces. This creates a far flow velocity v according to equation 
JIJ which, in the absence of any other particle, entrains particle 2 at r = (p, 0) with 
velocity U2 ~ v(r) given by equation (jfJJ). We now introduce particle 3 at r' = (p',0). 
It obstructs the flow and thus modifies u 2 . 

The particle being force-free, the lowest force-distribution moment it can exert is 
a force dipole, 

Sy(r') = A.faMo)- 1 ^ J a dv"{v" r') iVj (r"), (10) 

where the integration is over the surface of particle 3, and X s is a dimensionless factor 
of order 1. The force dipole, located at r', changes the flow velocity at r by 

<5<(r,r') = S jk (r')(a/wB )d J A kl (r - r'). (11) 

This holds regardless of confinement. In an unconfined system the far flows decay 
as 1/r, leading to S ~ (r')~ 2 and Sv s ~ (r')~ 2 |r — r'|~ 2 . Then, upon integration 
over the third-particle position r', one gets a correction ~ 1/r, which renormalizes 
the coefficient of the bare 1/r interaction. In Q2D, however, the situation is different. 
First, we realize that, due to the vanishing of v z and the symmetry of the mid-plane, 
only the terms i,j = x,y of Sij in equation (|10|) are non-zero. Then, since the far 
flows decay as l/p 2 , one gets Sv s ~ (p') _3 |p — p'|~ 3 which, upon integration over p' , 
yields a correction ~ l/p 4 . At large distances this is much smaller than the bare pair 
interaction (~ l/p 2 ) and will be neglected. 

Since the force dipole created by particle 3 has a negligible far-field effect in Q2D, 
we proceed to the third moment of the force distribution, 

T ijk (v') = X^a/wBo)- 1 ^ J dr"(r" - v'W - r')^(r"), (12) 

where At is another dimensionless prefactor depending only on the confinement ratio 
a/w. In principle, At (and X s ) could be calculated using the short-range hydrodynamics 
in a Q2D geometry ^HI EH)- yet we shall not pursue this technically complicated 
calculation here. The vanishing of v z and the mid-plane symmetry make all terms 
Tijk, for which one or three of the indices are z, vanish. In addition, the terms with 
(i,j, k) — (x, y) are negligible at large distances compared to those with two z indices, 
since they involve two additional powers of l/p'. We are thus left with T zzi , i — (x, y), 
yielding through equations (JIJ and l|12fl 

T zzl (p>) = h t a 2 (l A^p'Mj. (13) 

In the integration we have assumed again a parabolic vertical profile, H(z/w) = 
l—4z 2 /w 2 ^Hl- Note that within this assumption the similar leading terms in l/p' from 
all higher moments of the force distribution, involving higher z derivatives, vanish. The 
moment T(p') changes the flow velocity at p by 

5vl(r,r') =T ZZJ (p')(a/wB )d zz \ z=0 A Jt (p - p',0). (14) 

Combining equations ©, H13JI and l|14|l. we find the correction to the velocity of 
particle 2 due to particle 3, 

5u 2i = (a/iuB )CiAy(p - p')A jk (p')f lk (15) 

C - - A — f 1 - —\ (l - —\ 

1 ' 3u; 2 \ 3w 2 J \ 5w 2 J ' 
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where the reference to z — is hereafter omitted for brevity. 

Equation i|15fl yields the correction to u 2 given a certain position p' of particle 3. 
We now wish to average this correction over all possible p', 

{Sua) = J d 2 p'p(p,p')8u 2 {p lP '), (16) 

where p{p, p') is the probability density of finding particle 3 at p' given that particle 
2 is at p and particle 1 is at the origin. Employing the superposition approximation, 
we take this probability as 

P (p, p>) ~ -^g(p')g(\ P - p'\) * -^[1 + h(p') + h(\ P - p'\)}, (17) 

where g(p) is the pair correlation function of the Q2D suspension (normalized to 1 at 
p — > oo), and h(p) = g(p) — 1. We have assumed that the monolayer of particles is a 
disordered 2D liquid, and thus the pair correlation has no angular dependence. 

Using equations (|15|l - l|17|) . and specializing to the relation between (Su2 X ) and 
fix, we obtain the average corrections to the coupling mobilities, 

5Bt iT (p) = {a/w)B Q 5A L , T {p) 

6A L , T (p) = (tto 2 )- 1 ^ J d 2 p'[l + 2h{p')][A xx (p - p')A xx (p') + 

A xy (p-p')A xy (p')}, (18) 

where the longitudinal and transverse couplings are obtained, as in section[31 by taking 
p = px and p = py, respectively. [Note that equation l(T%|l . for h(p) = 0, has exactly 
the same form as the one used in based on a less detailed analysis.] 

A brief examination of equation Ijl8|l raises two expectations. First, the integrand 
scales as (p')~ 2 \p — p'\~ 2 which, upon integration over p', is expected to yield a 
correction ~ 1/p 2 , i.e., a renormalization of the prefactor of the bare coupling in 
accord with hydrodynamic screening. Second, since the fields decay slowly with 
distance, one expects features (oscillations) in the static correlation function h(p) to be 
smoothed by integration and thus be manifest only weakly in the dynamic coupling. 
As is shown below, due to the unique shape of the far flow in Q2D, both of these 
expectations turn out to be false. 

The convolution integral of equation l |18| l is worked out in the Appendix. The 
result is 



5A LjT (p) = C(j)w 2 



(19) 



where C = — 2(w 2 /a 2 )AgCi, and the results for the longitudinal and transverse 
couplings are identical. Equation i|19|) is the main result of the current work. It 
gives the leading three-body correction to the pair hydrodynamic interaction in Q2D. 
If static pair correlations are neglected, h(p) = 0, this correction vanishes, as was 
presented in |l()l 111) . Moreover, the correction for a nonzero h(p), although finite, 
is short-range: since the suspension is disordered, the static pair correlation h(p) 
decays exponentially beyond an equilibrium correlation length. Thus, the correction 
in equation l|19fl decreases with distance much faster than the 1/p 2 decay of the bare 
interaction. This remarkable disagreement with the usual notion of hydrodynamic 
screening (i.e., concentration-dependent change of the prefactor at large distances), 
known from unconfined systems, is fully confirmed by our experimental measurements 
as is seen in figure |2{e). 
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The first term in equation l|19|) represents the leading correction to equations 
and (jHJ as the inter-particle distance p decreases. It scales as h(p)/p 2 and, therefore, 
should be in phase with the oscillations of the static pair correlation. The second term 
is smaller, of order ah(p)/p 3 . Higher-order corrections, which have not been treated 
here, come from several sources. Corrections of 0(a 2 / ' p A ) arise from the integration in 
equation (|18Jl (see Appendix), as well as the force-dipole terms (5v s ) discussed above. 
Terms of 0(h 2 /p 2 ) come from the approximation employed in equation 117J1 . Finally, 
as p becomes comparable to w, short-distance hydrodynamics set in and the far flow 
, on which our entire analysis has been relying, should be corrected U^J. Thus, 
the expression formulated in equation 119fl is relevant to suspensions of sufficiently 
high concentration, such that the static pair correlation h(p) is significant at distances 
p > w where the far flow holds. 

To compare equation Ijl9(l with experiment we need the static pair correlation 
function for the various values of area fraction. These functions are directly measurable 
from snapshots of the Q2D suspensions such as the one shown in figure da). (For 
more details, see HOI-) The results are presented in figure 



o 



0123456789 

p/w 

Figure 3. Static pair correlation function as measured from suspension 
snapshots. Plotted is h(p) = g(p) — 1, where g(p) is the pair correlation function 
normalized to 1 at p — > oo. The inter-particle distance is scaled by w. Area 
fractions are, from the bottom up, <j> = 0.254 (red), 0.338 (orange), 0.547 (green), 
and 0.619 (blue). Curves are vertically shifted by 1 for clarity. 

The measured h(p) are subsequently used in equation (jT§)l to calculate the 
correction to the pair interaction. The results are shown in figure 0] With the 
coefficient A of the bare interaction already found, there is only one unknown 
coefficient, C, fitted as C ~ 0.85. Working back the definitions of the various 
coefficients introduced during the calculation, we find the coefficient of the third force 
moment for our experimental system (a/w ~ 0.45) to be At ~ 1.77. Comparing 
the fits in figure 0] with those in figure [3 we see that the introduction of the 
concentration-dependent correction leads to a good agreement between the theory 
and the measured longitudinal interaction not only at asymptotically large distances 
(as was presented in ^]) but also at intermediate distances, p > 2w. As 
anticipated, equation lj*H?)l is particularly successful for the high-concentration systems 
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p/w p/w 




Figure 4. Concentration effect on the longitudinal coupling diffusion coefficient. 
The diffusion coefficient is scaled by Dqcl/w and the inter-particle distance by 
to. Area fractions are <f> = 0.254 (a, red), 0.338 (b, orange), 0.547 (c, green), and 
0.619 (d, blue). Dashed lines are a fit to equations |I71. I19I with A = 0.36 (already 
fitted in figure[5J and C = 0.85 for all panels. Insets focus on the deviation of the 
coupling from its long-distance behaviour. The oscillations at high concentration 
are in phase with those of the measured static pair correlation (see figure 151. 



[figure 0fc,d)] , where oscillations related to h(p) are clearly observed (compare with 
figure |3J) . The discrepancies seen in figure ^ between the theory and experiment for 
p < 2w can be attributed to several factors, as discussed above, which have not been 
included in the current analysis. 

Our calculation has yielded an identical correction to the transverse pair 
interaction, <5At = SA^. This does not agree with the measured transverse 
interaction, which does not exhibit a wiggly behaviour similar to that of Al (see 
figure |2J). Thus, we cannot currently offer an accurate analysis of the concentration 
effect on the transverse interaction at intermediate and short distances. A possible 
reason for the difference between the longitudinal and transverse modes may lie in the 
higher sensitivity of the latter to short-range hydrodynamics. To demonstrate this 
difference we show in figure \5\ the deviations of the longitudinal and transverse bare 
interactions from their large- distance behaviour in the limit of very small particles 
(a/w <C 1), as obtained from an exact solution JSj. The departure of the transverse 
interaction from its asymptotic behaviour is found to be more pronounced than that 
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Figure 5. Deviation of the pair interaction from its large-distance behaviour in 
the limit of very small particles. The deviation of the transverse interaction (T, 
dashed line) is larger than that of the longitudinal one (L, solid line). The curves 
were calculated using the exact solution for a/w <C 1 given in [Tgl . The large- 
distance behaviour in this limit is given by Al,t ±\/{p/w) 2 , with A = 9/16. 



7. Discussion 

The correlated dynamics of particles in Q2D suspensions are very different from those 
in unconfined systems. Confinement affects the decay of the dynamic correlation with 
distance (1/r 2 instead of 1/r) as well as its sign (the correlation transverse to the 
line connecting the particles becomes negative). In the current work we have focused 
on the effect of concentration, i.e., the presence of additional particles, on the pair 
interaction. Unlike unconfined suspensions, where the large-distance coupling changes 
with concentration, in Q2D suspensions increasing the concentration has no effect at 
large inter-particle distances (no hydrodynamic screening). The qualitatively different 
behaviour of Q2D suspensions has been theoretically and experimentally corroborated. 
It stems from the dipolar shape of the far flow induced by particle motion. As argued 
in ^3 , the key property of this flow is that it is governed by the displacement 

of liquid mass rather than the diffusion of liquid momentum. 

Concentration-dependent effects on the pair correlation do set in in Q2D 
suspensions, yet only at intermediate and short distances. The range of these 
corrections is determined by the larger of two lengths: the range of the static pair 
correlation h(p) of the suspension and the confinement width w. In the former 
case, which is valid for sufficiently high concentration, the spatial dependence of the 
dynamic coupling reflects the features of the static pair correlation. In this case we 
have been able to provide a quantitative account for the concentration effect on the 
longitudinal interaction at intermediate distances, which is in good agreement with 
the experimental measurements. In the latter case (range of h shorter than w), the 
correction depends on short-range hydrodynamics whose analysis lies beyond the scope 
of the current work. 
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We have not treated the effect of particle motion perpendicular to the bounding 
surfaces. Such fluctuations must exist in practice, yet our results suggest that they 
have a minor effect at large and intermediate distances. This is expected since the 
flows produced by perpendicular fluctuations decay exponentially with distance and 
thus have only a short-range effect |18|. 
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Appendix 

We need to calculate the integral appearing in equation (|18|l : 

I(p) = J d 2 p'[l+2h(p')][A xx (p')A xx (p-p')+A xy (p')A xy (p-p%(20) 

where (omitting the prefactor) 

A xx (p) = (x 2 - y 2 )/p\ A xy (p) = 2xy/p 4 . (21) 
Changing to polar coordinates, p = (p, tp) and p' — (p' , tp'), we rewrite the integral as 

HP) = dp [1+2%) , , .o ZZJZi 3T2 A 22 ) 



p l2 [p 2 + p' 2 — 2pp' cos(tp' — ip)] 2 

from which it is evident that I(p) = I(p) is independent of tp. Hence, we may set 
(,5 = 0. This isotropy immediately implies also that the corrections to the longitudinal 
and transverse couplings, proportional to I(px) and I{py), respectively, are identical. 

The calculation of I(p) via equation l|22|) is a bit tricky for reasons similar to 
those encountered in electrostatics. We exclude two small areas ~ a 2 around the 
(integrable) singularities at p' — and p' = p, and divide the integration into three 
domains: (i) p' £ (a,p — a), tp' e (0, 2tt); (ii) p 1 £ (p + a, oo), (p' £ (0,27r); (hi) 
p' £ (p — a, p + a), tp' £ (a/p, 2tt — a/p). The contribution from the inner domain 

(i) vanishes upon integration over tp' for any p. (Note that the static pair correlation 
h(p) does not have an angular dependence.) Integration over tp' in the outer domain 

(ii) gives 27r[l + 2h(p')]/ (p') 3 . Subsequent integration over p' in domain (ii) yields 
2n dp'[l + 2h(p')]/(p') 3 + 0(a/p 3 ). However, due to the integrable singularity at 
p' = p, the intermediate narrow domain (iii) has a finite contribution as well, of 
opposite sign: — 7r[l + 2h(p)]/p 2 + 0(a/p 3 ). 

Thus, in the absence of static correlations, h(p) = 0, the leading terms in small 
a/p from domains (ii) and (iii) cancel and I{p^> a) =0. This result has already been 
obtained in A more detailed inspection reveals that the 0(a/p 3 ) corrections 

from domains (ii) and (iii) cancel as well and, hence, I — 0(a 2 / p A ). 

In the presence of static correlations, h ^ 0, we finally get 



I{ P ) = -27T 



hp) „r d m 



p 2 



0(a 2 lp% (23) 
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